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Abstract 

In the framework of the fey-factorization approach, we study the production of prompt 
photons and electroweak gauge bosons in high energy proton-(anti)proton collisions at mod- 
ern colliders. Our consideration is based on the amplitude for the production of a single 
photon or W /Z° boson associated with a quark pair in the fusion of two off-shell gluons. 
The quark component is taken into account separately using the quark-gluon scaterring 
and quark-antiquark annihilation QCD subprocesses. Special attention is put on the con- 
tributions from the quarks involved into the earlier steps of the evolution cascade. Using 
the Kimber-Martin-Ryskin formalism, we simulate this component and demonstrate that it 
plays an important role at both the Tevatron and LHC energies. Our theoretical results are 
compared with recent experimental data taken by the D0 and CDF collaborations at the 
Tevatron. 



PACS number(s): 12.38.-t, 12.38.Bx 



1 Introduction 

The theoretical and experimental studying the prompt photon and electroweak gauge 
boson production at high energies provide an important information about the nature of 
both the underlying electroweak interaction and the effects of Quantum Chromodynamics 
(QCD). In many respects these processes have become one of most important "standard 
candles" in experimental high energy physics [1—11] . 

In the previous publications [12-14], we have considered the production of prompt pho- 
tons and electroweak gauge bosons W ± and Z° in the fe t -factorization approach. Making 
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use of the /cy-factorization is motivated by the fact that it provides solid theoretical grounds 
for the effects of initial gluon radiation and intrinsic parton transverse momentum k?- We 
pay attention to individual contributions from the different partonic subprocesses. The idea 
of [12, 14] was in reexpressing the quark contributions in terms of gluon contributions, thus 
reducing the problem of poorly known and poorly calculable unintegrated quark densities 
to much better investigated gluon ones. Our studies, however, reveal the fact that the non- 
reducible quark distributions are of major importance for the processes under study. The 
goal of this paper is to clarify this point in more detail and to accomplish the calculations 
presented in [12, 14] by including the contributions which yet have not been taken into 
account. 

The outline of our paper is following. In Section 2 we recall shortly the basic formulas 
of the fey-factorization approach with a brief review of calculation steps. In Section 3 we 
present the numerical results of our calculations. The central point is discussing the role of 
each contribution to the cross sections. Section 4 contains our conclusions. 



2 Theoretical framework 

2.1 The subprocesses under consideration 

Our approach is the following. The starting point of consideration is the leading order 
0(a) and 0(aa s ) subprocesses: q + g* — > 7 + q, q + q 7 + g and q + q' — > W ± /Z°. 
These subprocesses are strongly depend on the unintegrated quark distributions in a proton 
f q (x, ky, /i 2 ). In contrast to the study [15] where the hard matrix elements of these subpro- 
cesses have been convoluted with the relevant unintegrated quark and/or gluon distributions 
in a proton, we try to reexpress the unintegrated quark densities in terms of gluon ones. Our 
main idea is connected with the separation of the unintegrated quark distributions into sev- 
eral parts which correspond to the interactions of valence quarks f^ v '(x, k 2 ^, /i 2 ), sea quarks 
appearing at the last step of the gluon evolution fjf'(x, k 2 -, /x 2 ) and sea quarks coming from 
the earlier gluon splittings fg S \x, k 2 -, fi 2 ) (see Fig. 1). In our approach, we simulate the last 
gluon splittings by the higher-order 0(aa 2 ) off-shell (i.e. /cy-dependent) matrix elements, 
namely g* + g* — > , y/W ± /Z° + q + q' . In this way we take into account the contributions from 
the f^(x, k 2 -, 11 2 ). To estimate the contributions from the f^(x, k 2 -, /i 2 ) and fi s \x, k 2 -, /1 2 ) 
we use the specific properties of the Kimber-Martin-Ryskin (KMR) scheme [16] which en- 
ables us to discriminate between the various components of the unintegrated quark densities^. 
Thus, the proposed scheme results to the following partonic subprocesses: 

g* + g* -> 1 /W ± /Z° + q + q', (1) 

q W + g *^ 1 / W ± /Z° + q\ (2) 

,W + g * _> 1 /W ± /Z° + q', (3) 

q + q^l + g, q + q' -> W ± /Z°. (4) 

To be precise, the gluon-gluon fusion subprocess (1) replaces the q^ +q^ annihilation, and 
the valence and sea quark-gluon scattering (2) and (3) replace the q( v > + q^ and q^ + q^ 

1 Below we will refer to the (x, , p?) contribution as to "reduced sea" component. 
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annihilation mechanisms. In the last two cases both valence and "reduced sea" q^ quark 
components are included. Of course, all incoming gluons in (1) — (4) are off-shell. To avoid 
the double counting we have not considered here q + q' — > W ± /Z° + g subprocess. 

As it was mentioned above, the proposed scheme was applied already to the prompt pho- 
ton [12] and electroweak boson [14] production at the Tevatron and LHC energies. However, 
in [12] we have neglected the "reduced sea" contribution. In [14], to estimate the "reduced 
sea" component only the contribution from the q^ + — > W ± /Z° subprocess has been 
taken into account since this contribution is the dominant one at the Tevatron energies (see 
also discussion in Section 3). In the present study we give a more accurate analysis of all 
possible contributions to the cross sections of processes under consideration^] and clarify the 
role of missing contributions to the "reduced sea" component at the LHC. 



2.2 Cross section for the inclusive r y/W ± /Z° production 

To calculate the cross section of the prompt photon and/or electroweak boson production 
in the framework of the /^-factorization approach one should convolute the off-shell matrix 
elements of subprocesses (1) — (4) with the relevant unintegrated quark and/or gluon dis- 
tributions. The contribution to the inclusive j/W^/Z production cross section from the 
off-shell gluon-gluon fusion (1) can be written as 

\M(g* + g* ->• 7 /W ± /Z° + q + g')' 2 



a(p + p -> 1 /W ± /Z° + X) = J2j 



kiTjA* )fg( x 2, k 2T , fx )dk 1T dk 2T dp 1T p 2T dydyidy 



2567r 3 (x 1 a; 2 .s) 2 

d(j>i d(j>2 dipi dip 2 



x 



(5) 



2tt 2tt 2tt 2tt 



where \M(g* + g* -> j/W ± /Z° + q + q')\ 2 is the off- mass shell matrix element squared (and 
averaged over the initial gluon polarizations and colors); yfs is the total energy of the process 
under consideration; krr, k2T, 4>i and 02 are the transverse momenta and azimuthal angles of 
the initial off-shell gluons (having the fractions X\ and X2 of the incoming protons longitudinal 
momenta); and y are the transverse momentum and rapidity of the produced prompt 
photon or vector boson; p 1T and p 2 T the transverse momenta of the co-produced quark and 
antiquark; y%, y 2 , and ip 2 are the quark rapidities and azimuthal angles, respectively. The 
formulas for the partonic subprocesses (2) — (4) are similar and can be written as follows: 

\M(q + g* ^^^/Z + q'Y 2 



■x 



16tt(xiX 2 s) 2 



x/ g (xi, k 2 T , fJ 2 )f g (x 2 , k^ T , fx 2 )dkl T dkl T dp^dydy' 



2vr 2vr 



\M(q + q^-f + g 



2 

-X 



167r(a;ia;2s) 2 



xf q (xi, k 2 1T , fi 2 )f q {x 2 , k 2 2T , ^ 2 )dk 2 1T dk 2 2T dp 2 T dydy g ^^ 



2 In the calculations of the prompt photon cross sections we will neglect the contributions from the so-called 
fragmentation mechanism [17]. It is because after applying the isolation cut (see [7-11]) these contributions 
amount only to about 10% of the visible cross section. The isolation requirement and additional conditions 
which preserve our calculations from divergences have been specially discussed in [12, 13]. 
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a(p + p->W ± /Z° + X)=Y / I 7 ^ , 2 \M(q + q' ^W ± /Z°)\ 2 x 

q J \XiX 2 S) 

x/ ff (xi, k 2 T) fi 2 )f q (x 2 , k 2 T , fi 2 )dk 2 1T dkl T dy^p-^, 

where y' and y s are the rapidities of the final quark or gluon. The analytic expressions 
for the off-shell matrix elements of subprocesses (1) — (4) has been derived in our previous 
papers [12, 14] (see also [18]). We only mention here that, in accord with the fc^-factorization 
prescription [19, 20], the off-shell gluon spin density matrix has been taken in the form 

£e^K^ = %^. (9) 

In all other respects our calculations follow the standard Feynman rules. If we average these 
expressions over <pi and 02 and take the limit k 2 T — > and k\ T — > 0, then we recover the 
relevant formulas in the leading order collinear approximation of QCD. 



2.3 The KMR unintegrated parton distributions 



In further analysis below we will use the unintegrated quark and gluon densities in a pro- 
ton which taken in the KMR form [16]. The KMR approach is the formalism to construct 
the unintegrated parton distributions f a {x, k 2 -, fi 2 ) from the known conventional parton dis- 
tributions xa(x, /j, 2 ), where a = g or a = q. In this approximation, the unintegrated quark 
and gluon distributions are given by [16] 



/ g (x,k|,^) = r g (k|,/, 2 )^ 



X 



x J dz [P qq (z)- z q g, 4) 9 (A - z) + P qa {z)- z g (~ k 



(10) 



fg(X, krp, /i ) — Tg(\irp, [1 



a s (k\) 



2tt 



dz 



(z)-q 

z 



X 



k^ 



+p t 



(*)-9 

z 



X 



k^ 



Q(A-z) 



where P a b{ z ) are the usual unregulated LO DGLAP splitting functions. The theta functions 
which appear in (10) and (11) imply the angular-ordering constraint A = ///(/i + |ky|) 
specifically to the last evolution step to regulate the soft gluon singularities. For other 
evolution steps, the strong ordering in transverse momentum within the DGLAP equations 



automatically ensures angular ordering^. The Sudakov form factors T 9 (ky, fjr) and T 9 (k^, fjr) 
which appears in (10) and (11) enable us to include logarithmic loop corrections to the 
calculated cross sections. The nonlogarithmic corrections can be taken into account by 
using the ^-factor [15, 23] K(q + q' ->• W/Z) ~ exp [C F vra s (/i 2 )/2] with C F = 4/3 and 
2 -p' /3 m 2 / 3 



3 Numerically, in (10) and (11) we have applied the recent Martin-Stirling-Thorne-Watt (MSTW) LO 
parametrizations [21] of the collinear parton densities a(x,/i 2 ). This choice is differs from the one [12-14] 
where the Gliick-Reya-Vogt (GRV) parton distributions [22] have been used. 
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The function f q (x, k^, /i 2 ) in (10) represents the total quark distribution function in a 
proton. Modifying (10) in such a way that only the first term is kept and the second term 
omitted, we switch the last gluon splitting off, thus excluding the fjf\x, k^, /i 2 ) component. 
Taking the difference between the quark and antiquark densities we extract the valence quark 
component ffj"\x, k 2 -, fi 2 ) = f q (x, k 2 -, fi 2 ) — f q (x, k 2 -, fi 2 ). Finally, keeping only sea quark in 
first term of (10) we remove the valence quarks from the evolution ladder. In this way only 
the fff\x, k|,, ii 2 ) contributions to the f q (x, k 2 ^, fi 2 ) are taken into account. 

The multidimensional integration in (5) — (8) has been performed by the means of Monte 
Carlo technique, using the routine Vegas [24]. The full C++ code is available from the 
authors on request^. This code is practically identical to that used in [12-14]. 

3 Numerical results 

We are now in a position to present our numerical results. In all our calculations according 
to (5) - - (8) the light quark masses were set to m u = 4.5 MeV, m d = 8.5 MeV and 
m s = 155 MeV. We have checked that the uncertainties coming from these quantities are 
negligible compared to the uncertainties connected with the scales in the unintegrated parton 
densities and strong coupling. We set m c = 1.4 GeV, mw = 80.403 GeV, mz = 91.1876 GeV, 
sin 2 6w = 0.23122 and use the LO formula for the strong coupling constant a s (fi 2 ) with 
rif = 4 active quark flavors at Aqcd = 200 MeV (so that a s (M§) = 0.1232). As it is often 
done, we choose the renormalization and factorization scales to be equal: hr = fip = yU = m*r, 
where is the transverse mass of the produced vector boson. In the case of prompt photon 
production we set the scale \i to be equal to the photon transverse energy BJj.. We will not 
study here the scale dependense of our results. This issue is addressed in our previous 
papers [12-14]. 

3.1 Role of the quark contributions 

We begin the discussion by presenting a comparison between the different contributions 
to the •y/W ± /Z° cross sections. In Fig. 2 we plot our results for the cross sections as 
a function of produced photon or gauge boson center-of-mass rapidity y. Here, we have 
performed the calculations for both the proton-antiproton and proton-proton interactions at 
the Tevatron and LHC energies, respectively. In the case of electroweak boson production, 
the cross sections are multiplied by the branching fractions f(W — > Iv) and f(Z — > l + l~). 
We set these branching fractions to f(W -> Iv) = 0.1075 and f(Z ->■ l + l~) = 0.03366 [25]. 
The additional cuts \y\ < 2.5 and \y\ < 4 have been applied in the case of prompt photon 
production at the Tevatron and LHC. The solid, dashed and dotted histograms in Fig. 2 
represent the contributions from the g* + g* — > r y/W ± /Z° + q + q' , q v + g* ■y/W ± /Z° + q' 
and q v + q' v — >■ W ± /Z° (or q v + q v — >• j + g) subprocesses, respectively. The thick solid 
histograms represent the sum of all contributions. One can see that the gluon-gluon fusion 
is an important photon production mechanism at both the Tevatron and LHC conditions. 

4 lipat ov ©theory, sinp.msu.ru 

5 In the present analysis the contributions from all 2 — s- 1 subprocesses have been corrected by the factor 
of 3/2 compared to those from [14]. We thank V.A. Saleev for drawing our attention to this point. 



5 



Moreover, it gives a main contribution to the cross section at the LHC. In the case of W ± /Z° 
production, the role of gluon-gluon fusion subprocess is negligible at the Tevatron and is 
increased greatly at the LHC energy: it contributes only about one or two percent to the 
total cross section at GeV and more than 40% at y/s = 14 TeV. The Compton-like 

subprocesses q v + g* — > r y/W ± /Z° + q' are also important. 

In Fig. 2, the dash-dotted histograms represent the "reduced sea" component. We find 
that this component gives approximately 30% contribution to the total cross section of 
prompt photon production at the Tevatron and approximately 20% contribution at the 
LHC. In the case of electroweak boson production, it contributes about 50% and 40%, 
respectively. As it was noted above, this component contains the + g* — > r )/W ± /Z° + q' , 
q (s) + g(s) _^ 7 + g ( Qr q (s) + q(s) _^ W ± /Z°) and gW + q^ -> 7 + g (or q^ + q^ W ± /Z°) 
subprocesses. The relative contributions are shown in Fig. 3. Note that thick solid histograms 
in Fig. 3 corresponds to the dash-dotted histograms in Fig. 2. Since all these subprocesses 
are mainly due to the quarks emerging from the earlier steps of the parton evolution rather 
than from the last gluon splitting, one can conclude that the quarks constitute an important 
component of the parton ladder, not negligible even at the LHC energies and not reducible 
to the gluon component. 

3.2 Comparison with the Tevatron data 

Now we turn to the comparison of our theoretical predictions with the experimental data 
on the prompt photon and W ± /Z° boson production at Tevatron. The data [7-11] on the 
inclusive prompt photon hadropro duct ion come from both the D0 and CDF collaborations. 
The D0 [7, 8] data were obtained in the central and forward pseudo-rapidity regions for 
two different center-of-mass energies, namely y/s = 630 GeV and y/s = 1800 GeV. The 
central pseudo-rapidity region is defined by the requirement |?7 7 | < 0.9, and the forward one 
is defined by 1.6 < |?? 7 | < 2.5. The more recent CDF data [9, 10] refer to the same central 
kinematical region |?7 7 | < 0.9 for both beam energies yfs = 630 GeV and y/s = 1800 GeV. 
The data on the electroweak gauge boson production are also come from the D0 [2-6] and 
CDF [1] collaborations. 

The results of our calculations are shown in Figs. 4 — 9. Fig. 4 confronts the double dif- 
ferential cross sections da / dE T drj of the prompt photon production calculated at yjs = 630 
and 1800 GeV in different kinematical regions with the D0 [7, 8] and CDF [9, 10] data. The 
solid histograms represent calculations in the scheme described above. For comparison, we 
also show (as dashed histograms) the predictions based on the simple 2 — > 2 QCD subpro- 
cesses with all quark components summed together. One can see that the both approaches 
give the very similar results which agree reasonably with the Tevatron data on the prompt 
photon cross sections within the experimental uncertainties. This fact demonstrates that 
the high-order corrections for prompt photon production connected with the terms not con- 
taining large logarithms are rather small. Our predictions based on the subprocesses (1) - 
(4) are rather similar to ones [26] based on the collinear QCD factorization with the NLO 
accuracy. 

Concerning the electroweak gauge boson production, the situation is slightly different. 
Figs. 5 and 7 display a comparison between the calculated differential cross sections da/dpr 
and the D0 and the CDF experimental data [1, 3, 4] at low p T (pr < 20 GeV) and in 
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the full pt range. These data have been obtained at y/s = 1800 GeV. In Figs. 6, 8 and 
9, we show the normalized differential cross section (1 / a) da / dpx and (1/cr) da/d\y\ of the 
W ± and Z° boson production in comparison with the data. The predictions based on the 
scheme (1) - - (4) are about a factor of 1.25 higher than the ones based on the simple 
2 — > 1 subprocesses. The main difference between the predictions is observed in the low 
Pt region. This difference can be attributed to the terms not containing large logarithms 
and connected with using of the high-order off-shell matrix elements mentioned above. In 
contrast with the prompt photon production, such terms are significant for the case of 
W ± /Z° production. In [15], an additional factor of about 1.2 was introduced ad hoc to 
eliminate the visible disagreement between the data and theory. The origin of this extra 
factor has explained [15] by the fact that the input parton densities (used to determine the 
unintegrated ones) should themselves be determined from data using the appropriate non- 
collinear formalism. The results of our calculations based on the scheme (1) — (4) show no 
need in this extra factor. Some overestimation of our predictions over experimental data at 
small pt can be, in principle, connected with problem of applicability of BFKL-like evolution 
in this region and requires an additional study. Note that the traditional QCD calculations 
(valid in a wide pt range) combine fixed-order perturbation theory (at present, up to NNLO 
terms [27-31]) with analytic soft-gluon resummation [32,33] and some matching criterion. 
These calculations gives a similar description of the data [1-6]. 

Additionally we have studied the effects of the non- logarithmic loop corrections to the 
gauge boson production amplitude. To do this, we have repeated the calculations based on 
the 2 — > 1 quark- ant iquark annihilation with the omitted i^-factor. The dotted histograms 
in Figs. 4 — 9 correspond to the results of these calculations. We have found a significant 
(by a factor of about 1.5) reduction of the predicted cross sections. Also, in our numerical 
calculations we have tested the rather old GRV (LO) parametrizations [22] of the collinear 
parton distributions in a proton (not shown in plots). We have found that the difference 
between the results based on the GRV and recent MSTW [21] collinear parton densities is 
negligible. 

Finally, we would like to mention that an additional possibility to distinguish the two 
calculation schemes comes from studying the ratio of the W ± and Z° boson cross sections. 
In fact, since W ± and Z° production properties are very similar, as the transverse momen- 
tum of the vector boson becomes smaller, the radiative corrections affecting the individual 
distributions and the cross sections of hard process are factorized and cancelled in this ratio. 
Therefore the results of calculation of this ratio in the scheme (1) — (4) (where the 0(aa s ) 
and 0(aa 2 s ) subprocesses are taken into account) and the predictions based on the 0(a) 
quark- ant iquark annihilation should differ from each other at moderate and high pt values. 
This issue have been studied in our previous paper [14]. 

4 Conclusions 

We have studied the production of prompt photon and electroweak gauge bosons in 
hadronic collisions at high energies in the /c r -factorization approach of QCD. The central part 
of our consideration is the off-shell gluon-gluon fusion subprocess g*+g* — > / ~y/W ± /Z + q+q'. 
The contribution from the quarks has been taken into account additionally. 
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To study the individual contributions from the different partonic subprocesses we have 
used the KMR scheme. We find that the gluon-gluon fusion is an important production 
mechanism of prompt photons at both the Tevatron and LHC conditions. At the LHC, 
it gives the main contribution to the cross section. In the case of W^/Z production, it 
contributes only about one or two percent to the total cross section at the Tevatron and 
more than 20% at the LHC energy. 

We demonstrate that an important contribution to the total cross sections of the processes 
under consideration also comes from the sea quark interactions. Notably, we found that the 
contribution of these subprocesses are mainly due to the quarks emerging from the earlier 
steps of the parton evolution rather than from the last gluon splitting. 
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Figure 1: Our approach to calculate the unintegrated quark distributions. 
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Figure 2: Differential cross sections of prompt photon or electroweak boson production at 
the Tevatron and LHC as a function of their cm. rapidity y. The solid, dashed and dotted 
histograms represent the contributions from the g* + g* — > r y/W ± /Z° + q + q', q v + g* — > 
•y/W ± /Z° + q' and q v + q v — > 7 + g (or q v + q v — > W ± /Z°) subprocesses, respectively. The 
dash-dotted histograms represent the "reduced sea" component. The thick solid histograms 
represent the sum of all contributions. 
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Figure 3: Different contributions from the "reduced sea" component to the prompt photon or 
electroweak boson production cross sections at the Tevatron and LHC. The solid, dashed and 
dotted histograms represent the contributions from the q s + g* — > r y/W ± /Z° + q', q s + q v — > 
7 + g (or q s + q' v — > W^ ± ) and q s + q s — > •y/W ± /Z° + g subprocesses, respectively. The thick 
solid histograms represent the sum of all contributions. 
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Figure 4: Double differential cross section da / dETdr] for the inclusive prompt photon 
hadroproduction calculated for \rj\ < 0.9 (a) and 1.6 < \rj\ < 2.5 (b) at y/s = 630 GeV 
and for \rj\ < 0.9 (c) and 1.6 < \rj\ < 2.5 (d) at = 1800 GeV. Solid histograms represent 
calculations in the "decomposition" scheme where all contributions described in the text are 
taken into account. Dashed histograms correspond to the predictions based on the simple 
2 — > 2 quark-gluon QCD interaction and quark-antiquark annihilation subprocess with all 
quark components summed together. The experimental data are from D0 [7, 8] and CDF [9, 
10]. 
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Figure 5: Transverse momentum distribution of the W ± boson production calculated at 
\fs = 1800 GeV. Solid histograms represent calculations in the "decomposition" scheme 
where all contributions described in the text are taken into account. Dashed histograms 
correspond to the predictions based on the simple 2 — > 1 quark-antiquark annihilation sub- 
process with all quark components summed together. Dotted histograms correspond to the 
simple 2 — > 1 quark-antiquark annihilation subprocess without .fT-factor. The cross sections 
time branching fraction f(W — > Iv) are shown. The experimental data are from D0 [4]. 




Figure 6: Normalized transverse momentum distribution of the W ± boson production cal- 
culated at y/s = 1800 GeV. Notation of the histograms is the same as in Fig. 5. The 
experimental data are from D0 [2]. 
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Figure 7: Transverse momentum distribution of the Z° boson production calculated at 
y/s = 1800 GeV. Notation of the histograms is the same as in Fig. 5. The cross sections 
time branching fraction f(Z—> l + l~) are shown. The experimental data are from D0 [3] 
and CDF [1]. 




Figure 8: Normalized transverse momentum distribution of the Z Q boson production cal- 
culated at \/s = 1960 GeV. Notation of the histograms is the same as in Fig. 5. The 
experimental data are from D0 [6]. 
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Figure 9: Normalized rapidity distribution of the Z° boson production calculated at y/s = 
1960 GeV. Notation of the histograms is the same as in Fig. 5. The experimental data are 
from D0 [5]. 
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